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We study the effects of random positional disorder in tlie transmission of waves in a ID Kronig- 
Penny model. For weak disorder we derive an analytical expression for the localization length and 
relate it to the transmission coefficient for finite samples. The obtained results describe very well 
the experimental frequency dependence of the transmission in a microwave realization of the model. 
Our results can be applied both to photonic crystals and semiconductor super lattices. 
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I. INTRODUCTION 

In recent years there is a high activity in the study of 
wave (electron) propagation through one-dimensional pe- 
riodic structures (see, for example, Ref. [l| and references 
therein). Much is already known about band structures 
of perfectly propagating waves in strictly periodic and 
relatively simple devices, and one of the current inter- 
ests is the influence of random imperfections that are 
commonly present in real experiments. These imperfec- 
tions are originated, for example, from the variations 
of the medium parameters such as the dielectric con- 
stant, magneticpermeability, barrier widths or heights 

The analysis of scattering properties of periodic-on- 
average (when periodic systems are slightly affected by a 
disorder) models with various kinds of disorder is mainly 
related to numerical methods. It is obvious that giving 
important results for specific models and parameters, the 
numerical approaches can not serve as a guide for the un- 
derstanding of generic properties caused by disorder. In 
this paper we try to fill this gap in the theory by the 
derivation of the localization length for the ID Kronig- 
Penney model, relating it to the properties of transmis- 
sion through a finite number of disordered barriers. 

Our analytical results are compared with the experi- 
mental data obtained for a single-mode microwave guide. 
We show that in spite of the standard restrictions of an- 
alytical results (restricted to infinite samples and weak 
disorder), comparison between theory and experiment is 
quite good. This fact is highly non-trivial since the ex- 
perimental data are strongly infiuenced by absorbtion in 
the waveguide walls; effect that is also not taken into 
account analytically. 

Our study is relevant to other types of ID stratified 



media, for example, to electron transport through ran- 
dom supcrlattices 10] (disordered arrays of semiconduc- 
tor quantum wells/barriers) or acoustic waves in random 
layered media [ill- Also, similar properties of the trans- 
mission are expected to occur in the ID quantum Kronig- 
Penney model (with alternating rectangular wells and 
barriers). 

In Sect. II the model is specified and the transfer ma- 
trix equations are derived. In Sect. Ill the experimental 
setup is briefly discussed and the numerical simulations 
for the transmission coefficient are compared with exper- 
imental results for the case of an array of 26 cells and 
different amounts of positional disorder. In Sect. IV we 
present the main experimental results and discuss some 
of the properties of transmission. In Sect. V we derive, 
for the regime of weak disorder, the analytical expression 
for the logarithm of the transmission in connection with 
the inverse localization length. We compare there too 
the numerical simulations and analytical results with the 
experimental data and show the effectiveness of our ana- 
lytical approach. In Sect. VI, we summarize our results. 



II. MODEL 

We consider an array formed by two alternating di- 
electric slabs with refractive indices tIq and rif, placed in 
an electromagnetic metallic-wall waveguide of constant 
width w and height h, see Fig. [T] For convenience, the 
layers with the refractive index Ua (rib) shall be referred 
as to the a-layers (6-layers). The lengths of the nth 
a- and ^-layers are denoted, respectively, by da{n) and 
db{n). The positional disorder in our model consists in 
randomly varying lengths of only one type of layer, say 
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FIG. 1: (color online) Kronig-Penney model with teflon bars 
of constant length db{n) — db{n + 1) = db, and air spacings 
da{n) defined by Eq. (f2fT|). 



the a-layer, such that 

da{n) = da + o-?7("-)7 {da{n)) = da, db{n) = 4. (2.1) 

Here a is r.m.s deviation of da{n) and its variance. 
Hence, r](n) is a sequence with zero average and unit 
variance. In this work we assume that r/(n) is random 
uncorrelated , i. e. 



{i]{n)ri{n')) 



ivH) = 0. 



(2.2) 



The angular brackets (. . .) stand for a statistical average 
over different realizations of randomly layered structure. 
Note that the random structure is periodic on average 
with the period d = da + db. 

In this work we shall treat the lowest TE mode of fre- 
quency v whose electric field E is defined by 



sm{nz/w)'^{x), 



0, (2.3) 



(see Fig.[T|). Within every a- or 6-layer, the function ^'(a;) 
obeys the ID Helmholtz equation, 



(2.4) 



dx^ 

with the wave numbers 
27r 



+ Kb ^a,6(x) = 0, 



{c/2wy 



(2.5) 



In the incoming-outgoing wave representation, the solu- 
tion of Eq. (|2.4p for the nth elementary cell reads 



(2.6) 



*Q„(a;) = e^p{ika[x - Xa{n)]} 

+A~ exp{-i/ca[a; - a^a ("•)]} 
inside the a„-layer, Xa{n) < x < Xb{n), and 
^b„{x) = B:^ exp {ikb[x - Xb{n)]} 

+B~ eyip{~ikb[x — Xb{n)]} 



(2.7) 



inside the 6„-layer, Xb{n) < x < Xa{n + 1). Here 
and are complex amplitudes of the forward/backward 
traveling wave, the coordinates Xa{n), Xb(n) denote the 
left-hand boundaries of the a„ and 6„ layers, respectively. 

With the use of the continuity conditions for the wave 
fimction '^a b(x) and its derivative at the boundaries x = 



Xb(n) and x = Xa{n + 1) one can obtain the transfer 
relation for the amplitudes A^_^_^ and A^ of two adjacent 
cells, 



A. 



n+1 



Qilin) Qi2(n) 
Q2i{n) Q22{n) 




(2.8) 



The transfer matrix Q{n) has the following elements, 
Qii{n) = [cos{kbdb) + ia+s\n{kbdb)]cyi'p[ikada{n)] 

- Q22("); (2.9a) 
'3i2(") = ia-s\n{kbdb)eyi-p[ikada{n)\ 

= Q2i{n). (2.9b) 

Here the asterisk stands for the complex conjugation, and 
we introduced the parameters a±. 



1 /fca _|_ fcfe 

2i \kb ka 



1. 



(2.10) 



The determinant of Q{n) is equal to unit, deiQ{n) = 1. 
Note that the transfer matrix Q{n) differs from cell to 
cell only in the phase factor exp[ifcaC?a(n)]. 

The transfer matrix equation for the array of N cells 
with or without positional disorder is 



^^+1 ) ^ qN 



At 
A 



(2.11) 



where. 



= QiN)Q{N - l)..Q(n)...g(2)g(l). (2.12) 

All matrices Q{n) {n — 1,2,..., N) have the same form 
(I2.9p . only differing in the values of da{n). In our fol- 
lowing numerical simulations and experimental set-up we 
have A'^^^ = 0. Thus the transmittancc of N cells is 
given by 

EE \A+^,/At\' = IQfir^ - (2.13) 
where is the scattering matrix element in the relation 



A + 



qN cN 

^11 ^12 

qN cN 

^21 ^22 



At 



A 



N+l 



(2.14) 



In the case of no disorder, rj{n) = , the length of a- 
layer does not depend on the cell number n, da{n) — da- 
Therefore, the unperturbed transfer matrix Q^"-* is de- 
scribed by Eq. (|2.9p with da{n) replaced by the constant 
length da- As is known (see, e.g., Ref. [l|), the transmis- 
sion through N identical cells is expressed in closed form 
as 



(0) 



N 



^12 sin(K(i) 



(2.15) 
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FIG. 2: (color online) The microwave guide. Left: overview 
of the brass waveguide of perimeter P = 234 cm, with height 
h — 1 cm and width w = 2 cm. Right: an enlargement of 
the waveguide showing the two carbon pieces (black) and one 
of the electric dipole antennas. The other lies on a brass lid 
(not shown) inserted near the other carbon piece. The white 
pieces are the teflon bars. 

where n is the Bloch wave number defined by the disper- 
sion relation 

cos(K(i) = co8{kada) cos{ki,dh) — Q!_|_ sm(kada) siii{kf,di,). 

(2.16) 

Expression (|2.15p indicates that the transmission is 

perfect {s[^''^ — 1) for all N where Q^""* = or when 
sm{NKd) I sin(K(i) = 0. The former occurs at the Fabry- 
Perot resonances fcfjdf, = m7r in the 6-layers. The latter 
produces — 1 Fabry-Perot oscillations in each spectral 
band, associated with the total system length Nd. We 
shall refer to the resonances k\,dh — imr as "teflon reso- 
nances" since in the experiment the 6-slabs are made of 
teflon, whereas the a-slabs are just air. 



III. EXPERIMENTAL SETUP 

In Fig. [2] we show the experimental setup of a mi- 
crowave ring guide of height h = 1 cm, width w — 2 
cm and perimeter P — 234 cm. The waveguide con- 
sists of iV = 26 cells, where the 6— layers are pieces of 
teflon of length db — 4.078 cm and the refractive in- 
dex rib — \/2.08. Two electric dipole antennas connect 
to a network analyzer. Also shown are the two carbon 
pieces used to absorb the electric field at both ends of the 
waveguide and thus mimic two infinite leads connected 
to each side of the array of teflon pieces and air segments. 
The frequency range is 7.5 to 15 GHz, corresponding to 
wave lengths from 4 to 2 cm. This arrangement has been 
used to study the transport effects of single impurities 
in the photonic Kronig-Penny model [T^. Earlier, in an 
analogous model with metallic screws instead of teflon 
pieces, the microwave realization of the Hoftstadter but- 



terfly was studied. The same configuration (metallic 
screws) was used to investigate transport properties of 
on-site correlated disorder [ij) IHj [3 ■ 

Obviously this waveguide is not rectilinear; in the ex- 
periment the teflon pieces and air segments are not per- 
fect parallelepipeds: one side is longer than the other by 
5 percent. However since the perimeter (234 cm) is much 
larger than the wave length of the electric field even in 
the regime of the first mode, it is expected that the recti- 
linear model is a good approximation. In fact, as shown 
in detail in p^ . a good quantitative agreement is found 
by defining an effective length of the teflon pieces and 
the air segments that are 1.95 percent larger than the 
smaller side of the teflon pieces and air segments. For 
example, the length of the smaller (inner) length of the 
teflon pieces used in the experiments reported here is 4 
cm, so the effective width da we use in our calculations is 
4.078 cm. We remark that this value is found by best fit- 
ting and is the only fitting, good for the whole frequency 
range of all our results presented here. 

In Fig. [31 curve (a), we plot the experimentally mea- 
sured value of |S'f^| (in what follows, the transmission 
spectrum) for the 26-cells periodic array. Curve (c) is the 
transmission spectrum |>5'j2^^|, calculated according to 
Eqs dmni), (Pl^ and unperturbed Eq. Note that 

the experimental transmission spectrum is about 1/5 of 
the theoretical one. This decrease of the signal is due 
to absorbtion by the metal walls of the waveguide. The 
value of the transmission spectrum is roughly constant 
over the frequency range, and can be taken into account 
phenomenologically by introducing an absorbtion factor. 
However, it is not the purpose of our work here to study 
the absorbtion effects. Our question is the global fre- 
quency dependence of the transmission spectrum on the 
frequency, giving us a possibility to reveal resonance ef- 
fects and the role of disorder. Note that the band struc- 
ture of the spectrum remains practically the same in spite 
of a strong absorbtion (see discussion in Ref. [l^. ITsl. Ila| ) . 

Inspection of the transmission spectrum |5j°'^|ofthe 
perfectly periodic array. Fig. [3l curve (c), demonstrates 
the effect of the teflon resonances on the transmission 
bands. We see two types of bands; namely, bands 1,2,4, 
and 6 show the — 1 = 25 oscillations mentioned above, 
whereas bands 3, 5 and 7 are flat i\s[l^'^\ « 1) around 

the resonance since |Qi2^| is zero at the resonance, and 
very small in some neighborhood around it. For future 
reference we shall refer to the second type of bands as 
the "resonance bands". Clearly, these bands disappear 
when db 0, turning into those with A'^ — 1 oscillations 
occurring for delta-function potentials. 

In accordance with Eq. (|2.5p . the first (lowest) mode 
in the air spacings (a-slabs) opens at the cut-off fre- 
quency i^™* — {c/2wna) = 7.5 GHz while in the teflon 
b-layers it opens at z^J"* = {c/2wnb) = 5.2 GHz, which is 
less than i^™* {ua < rib). On the other hand. Fig. [3] 
specifies the bottom of the first transmission band at 
= 7.387 GHz. So that the real cutoff Per- 
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FIG. 3: (color online) Transmission spectrum for a 26-ceIls 
array. Top: 7.3 < f < 11 GHz. Bottom: 11 < i/ < 15 
GHz. Curve (a): experimental data for the periodic array. 
Curve (b): |5i2|, multiplied by 1/2, for a slightly random 
(e = 0.049) array. Curve (c): IS'Ja"'^! for the perfectly periodic 
array {da = df, = 4.078 cm). Dashed vertical lines mark the 
position of the teflon resonances. 



tains to the interval i/™* < < i-™* where the 

wave number fc^ in the air spacings is purely imaginary, 
ka — i(27r/c)-\/ (c/2u')2 — ly"^. The transmission in this 
regime is due to tunneling through the air spacings. This 
fact makes the profile of the first band (see Fig. [3l curve 
(c)) somewhat special: there is a dip in the transmission 
band right at the frequency where the first mode opens 
in the air. 

Note that except for the first band, the positions and 
widths of the transmission gaps as well as the band pro- 
files of the experimental curve, are weU reproduced by 
the theoretical calculations. The discrepancy for the first 
band is understood since for low frequencies the wave 
length is not sufficiently small compared to the perime- 
ter of the circular waveguide, and hence the rectilinear 
waveguide model fails. For the second and higher bands, 
the agreement is better and the experimental curve does 
show some evidence of the small iV — 1 band oscillations 
predicted by the model. Clearly, these do not appear as 
perfectly regular oscillations and this irregularity may be 
caused by experimental imperfections due to variations 
in the length and positions of the teflon pieces. 

The maximum deviation \di,{n) — dii\max in the length 
of the teflon pieces is about 0.01cm and the maximum 
deviation \da[n) — da\max in the air spacings, due to the 
placement of the teflon pieces, is estimated to be about 
0.04 cm. Are imprecisions of this order sufficient to break 
the regularity of the oscillations? 

To check this, we performed a simulation with a disor- 
dered sequence assuming \da{rL) — da\max — 0.04cm and 
an error described by a random sequence 77(71) with a uni- 
form random distribution in accordance with Eqs. (j2.ip . 
(1221), (Pr^ and 11113]). The result is plotted in Fig. H 
curve (b). This curve has been multiplied by 0.5 in order 
to show it together with the perfectly regular case and 



case 


e/10-2 




(f) 


very weak 


0.49 


0.28 ■ 10-2 


8.0- 10-^ 


weak 


3.00 


1.77- 10-2 


3.0- 10-" 


medium 


12.30 


7.07- 10-2 


5.0- 10"^ 


strong 


49.00 


28.30 ■ 10-2 


8.0- 10-2 



TABLE 1: Parameter values of random disorder. Here d — 
da + db — 2da ~ 8.16 cm 

the experimental data. Inspection shows that indeed the 
assumed small error is enough to break the regularity of 
the band oscillations giving a better agreement with the 
experimental data of the supposedly regular array. 

IV. DISORDERED ARRAY 

Let us now move to intentionally disordered arrays, 
with the lengths of all teflon bars constant, db{n) = 
db while the air layers have random lengths given by 
Eqs. (|2.ip . (|2.2p . In our experimental and numerical cal- 
culations, the sequence 77(71) is an uncorrelated random 
function uniformly distributed in the interval [— a/S, VS], 
with unit variance. 

Apriori it is not known how large a random deviation 
from the average value da should be to observe weak, 
medium or strong disorder effects in the transmission. 
We tentatively classify the amount of disorder by the 
value of the maximum deviation from the average length 
of the air spacing divided by the average length of the 
cells, 

d 

Table H] shows the values of e we consider in this work, 
together with corresponding values of relative r.m.s. a/d 
and (a/d)'^. The latter quantity is needed to ease the 
comparison with analytical results obtained below. The 
case of e = 0.49 x 10^^ was discussed above to simulate 
the errors in the experimental setup. We call it the case 
of extremely weak disorder. Similarly, e = 3.0 x 10"^, 
12.3 X 10^^ and 49.0 x 10"^ cm, respectively, are called 
the weak, medium, and strong disordered cases. 

Fig. m shows the transmission for the array of 26-cells 
with the positional disorder. Compared with the case of 
weak disorder (Fig. [4^), for medium disorder (Fig. [JJd) 
only the first two gaps are clearly distinguishable; the 
third only partially. There is no trace of the TV — 1 
oscillations in the transmission bands, and the second, 
fourth, and sixth transmission bands have decayed sub- 
stantially. However, remnants of the resonance bands are 
still recognized, and so this can be considered the regime 
of medium disorder. For strong disorder (Fig. Uj:) the 
first two transmission bands have disappeared. There is 
no longer any evidence of the band structure of the un- 
perturbed array. But still the transmission spectrum is 
close to one in the vicinity of the teflon resonances. 
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FIG. 4: (color online) Transmission spectrum \Si2 \ for 26-cells arrays: (a) weak (e = 3.0 x 10~^), (b) medium (e = 12.3 x 10"^), 
and (c) strong disorder (e = 49.0 x 10"^). Transfer matrix calculations are shown by thin black solid curves, and experimental 
measurements by thick red solid curves. Dashed curves correspond to the analytical expression (|5.18p for the inverse localization 
length. Perpendicular lines mark the position of the teflon resonances. 



The results so far discussed pertam to the array of 
TV = 26 cells so the question arises about the effects on 
larger arrays. Our experimental setup does not allow 
for the implementation of much larger arrays. However, 
given that the numerical simulations are in good corre- 
spondence with the experimental data for = 26 cells, 
we now consider only numerically, and later analytically, 
larger arrays for the same three cases: weak, medium, 
and strong disorder. In Figs. - [St we plot for 
an array with TV — 100 cells, for weak (Fig. [5^), medium 
(Fig. (SJ)), and strong (Fig. [St) disorder. Similarly, in 
Figs. 5d - 5f we plot |S'f^| for an array with TV = 400 cells, 
again for weak (Fig 5d), medium (Fig. [St), and strong 
(Fig. [Sf) disorder. Comparing, for example. Figs. 5a and 
5d, corresponding to weak disorder for TV = 100 and 
TV = 400 cells, respectively, we see that the effect of in- 
creasing the size of the array, keeping the same amount 
of disorder, is to further decrease the transmission, con- 
sistent with the localization theory. However, this de- 
crease occurs only away from the teflon resonance fre- 
quencies. For medium and strong disorder, transmission 
has decayed bellow 10"'^ for most of the frequencies ex- 
cept around the teflon resonances. Thus the localization 
is not homogeneous at all: the teflon Fabry-Perot reso- 
nances strongly suppress the localization. 



V. LOCALIZATION LENGTH 

In this section we derive an analytical expression for 
the localization length Lioc and relate it with the exper- 
imental data. In the case of weak positional disorder. 



(5.1) 



an analytical expression for this quantity can be obtained 
as follows. First, we expand the transfer matrix Q{n) 
defined by Eq. (|2.9p . up to quadratic terms in the per- 
turbation parameter ka<Jrj{n)^ 



Q^n) « |l - (^ii^Ml^l Q(0) + fc„^^(„)Q(i). (5.2) 



The unperturbed Q*-"-* and first-order Q^^^ matrices are 
suitable to be presented in the form 



Q 



(0) 



Q 



(1) 



; (5.3a) 
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FIG. 5: (color online) Transmission spectrum j5'i2| (solid curves) and exp{~lioc) (dashed curves) for A'^ — 100 and TV — 400: 
(a) to (c) N=100 cells; a) weak, b) medium and c) strong disorder, d) to f) N=400 cells; d) weak, e) medium and f) strong 
disorder. 



u = [cos(fcf,(if,) + ia^ sin(fcb'^f))] exp{ikada) , (5.3b) 

V — ia^ sin(fcbC?6) <ixp{ikada) , (5.3c) 

det Q(°) = det Q(i) = \u\^ - \v\^ = 1. (5.3d) 

Also, it is useful for further calculations to take into ac- 
count that the real and imaginary parts, Ur = Rew and 
Ui = hnu, of the matrix element u can be expressed as 



Ur — cos(K(i), uf = sin^(K(i) + \v\'^ 



(5.4) 



The first equality is identical to the dispersion relation 
(|2.16[) . while the second one is a direct consequence of 
the matrix unimodularity (j5.3d|) . 

In order to extract the effects that are solely due to 
disorder, it is conventional to perform a canonical trans- 
formation to the Bloch normal-mode representation in 
the transfer relation 



A + 



4" 

A- 



P 



(5.5a) 
(5.5b) 



The transformation matrix P is specified in such a man- 



ner to make the unperturbed matrix Q^*^) diagonal. 



exp(-l-iKd) 







exp(— ZKO?) 



(5.6) 



in complete accordance with the Floquet theorem Il7j, 
or the same, the Bloch condition [l^. The solution of 
the problem for the eigenvectors and eigenvalues of Q^"^ 
results in 



P = 



i\/l3+ -iv*/f}. 



iv/ (3^ 



= 2 sin(K(i) [iti ^ sin(K(i)], 
f3l(3^_ = A\vfsm^{Kd), 

detP = detP-^ = \vf{P^^-PZ^) 



(5.7a) 

(5.7b) 
(5.7c) 
1. (5.7d) 



After substituting Eqs. (|5.3ap and (I5.7ap into 

the canonical transfer relation (|5.5ap , one can obtain the 
explicit perturbative recursion relations for the new com- 
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plex amplitudes, 



A + 



k^fj rf{n) ika(Trj{n)ui 



sin(K(i) 



A. 



n+l 



kacrr]{n)v* 
sin(K(i) 

k1(T'^rf{n) ika(J-q(n)ui 



1 - 



ka(jri{n)v 
sin(K(i) 



sin(K(i) 

exp(— 



exp(iK(i)A^ 

(5.8) 

exp(— 

(5.9) 



Now one can see from these equations that one equa- 
tion can be directly obtained from the other just by com- 
plex conjugation, if we suppose that — A^* . In other 
words, it is convenient to seek the amplitudes A^ in terms 
of action-angle variables. 



A^ = i?„exp(±i6i„) 



(5.10) 



In order to derive the equation for the real amplitude i?„, 
we muhiply Eq. ([5^ by Eq. ([SJ]! . Within the second 
order of approximation in the perturbation parameter 
kaorjin), we have 



R2 



1 



2kaari{n)\v\ 



' sin(K(i) 
kl(j'^rf{n) 



sin (kcT) 



uf + + 2u,\v\ cos(26l„ + kada)]. (5.11) 



The logarithm of Eq. (|5.1ip . that determines the local- 
ization length, is also be expanded within the quadratic 
approximation. 



In 



2kaari{n)\v\ 



??2 

i?^ J sin(K(i) 
2kla^r,^{n)\v\^ 



sin(20„ 



sin (nd) 



1 - sm\2e^ + kada) 



+ — cos(26'„ -I- kada) 
\v\ 



(5.12) 



Now we are in a position to write down the expression 
for the inverse localization length L^^^ that is known to 
be defined as follows flQll. 



1 



-"•n+l 



(5.13) 



The average (ab) is performed over the disorder ri(n) and 
the average ab is carried out over the rapid random phase 
9n- Within the accepted approximation and for uncor- 
related disorder, see Eq. (|2.2p . we may regard the ran- 
dom quantities 77(77.) and rf'{n) to be uncorrelated with 
trigonometrical functions, containing the angle variable 



On- Moreover, it can be shown (see, e.g., Ref. [i0|) that 
the distribution of phase within the first order of 
approximation in a weak disorder, is homogeneous (the 
corresponding distribution function is constant). There- 
fore, after averaging over 6'„ of Eq. (|5.12p . the term lin- 
ear in 77(77) and the last term in the brackets vanish and 
siv? [20 n + kada) is replaced with 1/2. As a result, we get 



, sin^(fcf,dfc) 
2d sin^ {nd) 



(5.14) 



This expression is in complete correspondence with that 
obtained in Refs. JJJ and 20] using a different approach, 
and reduces to Eq. (13) of Ref. [2l| for the limiting case of 
delta-like barriers. The appearance of the term s\T?{k},dh) 
in the numerator of Eq. (|5.14p indicates that at frequen- 
cies obeying the teflon resonance condition fcbdf, = ttztt, 
the localization length turns into infinity. That is, the 
random array becomes transparent and this is what is 
observed in the experimental and numerical transmission 
coefficient plotted in Figs. 3-5. 

As is known, the localization length is directly related 
to the transmittance T^r for a finite array of the length 
L = Nd, according to the famous relation (InTjv) = 
—2L/Lioc- In view of this relation and recalling that 
Tat — see Eq. (|2.13|) . it is convenient to introduce 

the resettled inverse localization length as 



{ln\S'^2\)^-L/Lio 



loc ' 



According to Eq. (|5.14p . one can get 

liol = i<^/d)'Fi,.)N. 
Here we introduced the form-factor 

sin^(fcbdb 



F{l^) = (kad)' 



2 sin (nd) 



(5.15) 



(5.16) 



(5.17) 



that specifies the frequency profile of the inverse localiza- 
tion length and is determined only by the parameters of 
the underlying regular array. The rescaled inverse local- 
ization length. (j5.16p . increases linearly with the number 
of cells N and quadratically with the amount of disorder 
a/d. 

In Fig. [5^, we plot the form factor F{i') for the pa- 
rameters of our system, da — db = 4.078 cm, tt.^ = 1, 
Ub = V2.O8. In Fig. [6b, we present the frequency de- 
pendence of for an array of 400 cells for the cases of 
weak, medium, and strong disorder. These figures show 
that the attenuation of transmission is larger for the 6th 
band and least for the first band, and that the states are 
completely extended at the teflon resonances. 

In order to compare our analytical results with the ex- 
perimental and numerical ones, we define the theoretical 
value IS'f^lth of the transmission spectrum as follows, 

IS^'.ltH ^ eM-liol) - cxp[-(a/d)2F(7.)A^]. (5.18) 

This is what is plotted in green dashed lines in Fig. 0] 
for the cases of weak, medium and strong disorder in 
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FIG. 6: (color online) (a) Form-factor F{v,da, di,,na,nt) for 
da — db = 4.078 cm, Ua = 1 and nt = V2.08; (b) inverse 
localization length for A'^ = 400 for weak (solid curve), 
medium (dashed curve), and strong (dotted curve) disorder. 



the = 26 cells array and in Fig. [5] for the array with 
N = 100 and N = 400. 

Inspection of Figs. 4-5 reveals that the theoretical ex- 
pression for the transmission spectrum provides a very 
good description for the case of weak disorder for the 
whole range of frequencies. For medium disorder, the 
agrement is moderate but only up to the region of the 
first resonance band. The higher the frequency the higher 
the discrepancy, in accordance with the range of validity 
(|5.ip of the expression (|5.16p . namely, (fcaa)^ <C 1. 

We would like to stress that the statistically aver- 
aged analytical results of this section required the struc- 
ture to be sufficiently long. Therefore, a better corre- 
spondence with numerics and experiment is expected for 
longer samples. The data in Figs. [3] and [H] confirm this 



fact. Indeed, in contrast with = 26 (Fig. |4]), for 
N = 100 and N = 400 (Fig. [5]) we can see much better 
agreement between theory and experiments/numerics for 
longer chains. The quantity exp(— Zj~J) gives a good fit 
to the numerical data even for the case of medium disor- 
der and up to J/ = 11.5 GHz. Observe also that already 
for N — 100 the non-resonance bands have completely 
disappeared for the medium disorder case. 



VI. SUMMARY 

We studied the effects of the uncorrelated positional 
disorder in the ID Kronig-Penney model paying main at- 
tention to the transmission for weak, medium, and strong 
disorder. We compare experimental data obtained in the 
microwave Kronig-Penney setup with the direct calcula- 
tions based on the transfer matrix method. We show that 
in spite of a very strong absorbtion in the metallic walls 
of waveguide, the experimental frequency dependence of 
the transmission is similar to that obtained numerically 
by neglecting the absorbtion effect. In particulary, the 
position of the TV — 1 Fabry-Perot resonances, as well as 
the teflon resonances are in correspondence with numeri- 
cal data. The structured similarity between experimental 
and numerical data can be explained by the different na- 
ture of the absorbtion in comparison with the resonance 
effects. Namely, the latter are due coherent effects, in 
contrast with the non-coherent nature of the absorbtion. 

Another part of our study is related to the expres- 
sion for the localization length obtained in the pertur- 
bative approach for an infinite chain of scattering barri- 
ers. Although in the experiment we have a relatively 
small number of 26 teflon barriers, we show that the 
expression for the localization length can be effectively 
used for the description of the transmission coefficient 
through the finite sample of barriers. To do this, we 
have used the relation between the transmission coeffi- 
cient and localization length that involves the finite size 
of samples, and applied it to the experimental situation. 
As a result, we have found that the analytical expression 
for the transmission coefficient reproduces quite well the 
frequency dependence of experimentally obtained data. 
Note again, that the global correspondence between ana- 
lytical and experimental results occurs in the presence of 
strong absorbtion that is non-avoidable experimentally. 
The observed correspondence between experimental and 
analytical results indicate that the expression for the lo- 
calization length is a working quantity for finite samples, 
a fact that is not commonly used in applications. 

The method we used to derive the expression for the 
localization length can be generalized to a more general 
case of the correlated disorder [l4, 15, 20, 21, 22, 23, 2|, 
[Hj HEI , for which anomalous effects in the transmission 
are predicted and experimentally observed. Our results 
can be applied to photonic crystals and electron super- 
lattices, as well as to propagation of acoustic waves in 
disordered systems. 
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